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Annual  progress  report  (2011-2012)  for  Contract  W81XWH-09-1-0726 
SYSTEMS  BIOLOGY  OF  GLUCOCORTICOIDS  IN  MUSCLE  DISEASE 


Introduction 

Duchenne  muscular  dystrophy  (DMD)  is  the  most  common  and  incurable  muscular  dystrophy  of  childhood. 
Muscle  regeneration  fails  with  advancing  age,  leading  to  considerable  fibrosis.  Corticosteroid  therapy  in  DMD 
has  been  shown  to  improve  muscle  strength  and  function  both  acutely  and  over  long-time  frames.  The 
molecular  mechanisms  explaining  the  beneficial  effect  of  glucocorticoids  is  unknown;  and  the  beneficial 
response  of  DMD  patients  is  enigmatic  due  to  the  deleterious  wasting  effects  chronic  glucocorticoids  have  on 
normal  muscle.  The  goal  of  this  proposal  is  to  use  computational  systems  biology  approach  to  model  the 
molecular  responses,  and  to  extend  and  validate  the  integration  of  signaling  and  transcriptional  networks  using 
our  pre-existing  high  throughput  mouse  and  rat  data  sets  of  transcriptional,  proteomic,  and  PK/PD  response  to 
glucocorticoids  and  human  neuromuscular  diseases.  We  also  use  laboratory  approaches  to  test  the  hypothesis 
that  the  effect  of  glucocorticoids  on  reducing  inappropriate  cross-talk  in  TGFp  networks  improves  muscle 
regeneration.  We  anticipate  that  the  computational  models  and  experimental  results  resulted  from  the 
proposed  research  will  provide  a  novel  model  for  efficacy  of  steroids  in  DMD. 


Body 

Progress  of  Specific  Aim  1 

Specific  Aim  1.  Provide  a  multi-scale  computational  model  on  the  integrated  acute  proteome  (signaling)  and 
transcriptional  response  regarding  metabolic  remodeling  and  TGFp  cascades  in  muscle  induced  by 
glucocorticoids  under  bolus  administration  of  glucocorticoids.  (August  2009  -  December  2010) 


Re-profiling  of  rat  acute  transcriptional  time  series  of  bolus  administration  of  glucocorticoids  using 
lllumina  BeadArray  gene  expression  chip.  (Year  2) 


In  order  to  effectively  infer  transcriptional 
networks,  we  have  re-profiled  the  samples  of  the 
rat  acute  transcriptional  time  series  of  bolus 
administration  of  glucocorticoids  using  lllumina 
gene  expression  BeadArray  chips.  The  cost  was 
covered  by  a  different  funding  resource.  The  old 
time  series  data  were  generated  using 
Affymetrix  U34A  chip  a  number  of  years  ago, 
this  outdated  chip  contains  only  5138  unique 
gene  (8740  probe  sets)  while  the  new  lllumina 
BeadChip  contains  21675  unique  genes/probes. 

From  our  initial  analysis  of  the  new  lllumina  time 
series,  we  have  found  many  genes  that  are 
affected  by  glucocorticoids  are  in  the  lllumina 
chip  but  not  in  the  old  Affymetrix  U34A  chip 
(Figure  1  shows  some  examples).  Although  both 
state  space  model  and  Bayesian  networks  we  used  in  our  network  reconstruction  approach  can  use  hidden 
variables  to  represent  possible  factors  contributing  to  regulatory  interactions  that  are  not  measured,  introducing 
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Figure  1.  Heat  map  of  some  transcripts  that  are 
contained  in  lllumina  gene  expression  chip  but  not  in 
Affymetrix  chip  and  are  significantly  affected  by 
glucocorticoids.  Red  color  indicates  up-regulated,  and 
green  for  down-regulated. 
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excessive  hidden  variables  may  significantly  increase  the  risk  of  model  over-fitting.  Given  a  larger  number  of 
measured  transcripts  in  the  new  time  series,  we  expect  that  we  will  obtain  more  reliable  and  complete  network 
results. 


Development  of  proposed  linear  state  space  model  for  genetic  network  reconstruction  (Year  1  and  2) 


We  developed  a  prototype  of  the  proposed  wavelet/linear  state  space  model.  As  described  in  the  proposal,  this 
approach  consists  of  two  components.  The  first  component  is  gene  modulization  by  wavelet  feature  extraction, 
extracting  transient  features  using  wavelet  transform  1-3  from  time  series  signals  and  modulizing  the  genes. 

The  second  component  is  molecular  network  construction,  in  which  we  estimate  the  structure  and  parameters 
of  a  molecular  network  to  quantify  the  interaction  relationships  among  the  identified  genes/gene  modules.  In 
the  proposal,  we  planned  to  use  two  alternative  approaches  for  the  second  component.  Approach  A  uses 
dynamic  Bayesian  network  (DBN)  to  estimate  the  parameters  and  conditional  probability  distributions  in  Eq.  (1) 
4-9  (reported  in  the  following  section).  Our  newly  proposed  Approach  B  uniquely  combines  wavelet  transform 
(trend  transform)  and  linear  regression  for  estimating  the  parameters  of  the  linear  regulation  model  (Eq.  (2)). 


Y,+1  =  GY,  +KU,  +  W,,(f  =  l,2,...r  -l) 


X,+1  =  GXf 


+  KU„Y,  =CX, +W,,(f  =  l,2,...r-1) 


(1) 

(2) 


Figure  2.  Simulation  experiment  on  genetic  network  reconstruction  using  Approach  B  wavelet/linear 
regulation  model,  (a)  Simulated  gene  expression  of  six  genes  and  an  external  stimulus,  (b)  The  regulatory 
relationship  among  the  genes  and  the  external  stimulus.  The  weight  of  the  lines  indicates  the  strength  of  regulation. 


We  developed  and  tested  the  combined  wavelet  analysis  and  linear  regression  method  (Approach  B)  on 
simulated  data,  and  obtained  promising  results.  We  generated  an  interaction  model  (Figure  2b)  consisting  of 
six  genes  (each  with  a  times  sequence  representing  its  change  in  response  to  an  external  stimulus)  according 
to  the  linear  Markov  model  with  given  parameters  and  an  external  stimulus.  We  simulated  different  types  of  the 
regulatory  relationship  among  the  genes  and  the  external  stimulus,  including  positive  stimulation  and  negative 
inhibition,  self-stimulation  and  self-inhibition,  and  gene  with  no  interaction  with  all  others  (Figure  2b).  The 
thickness  of  the  lines  indicates  the  regulatory  strength  (regulatory  coefficient).  Random  noises  were  also  added 
to  each  signal.  Then,  we  used  this  method  to  estimate  the  regulatory  coefficients  using  the  noisy  signals  (i.e. , 

G,  and  K  in  Eq.  (2)  for  wavelet/linear  regression).  Finally,  we  used  the  estimated  coefficients  to  reconstruct  the 
signals,  and  compared  them  to  the  original  signals.  The  closer  the  reconstructed  signals  (smooth  curves, 
without  noise)  are  to  the  original  signals  (fluctuating  curves,  with  noise  embedded),  the  better  the  parameters 
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are  estimated.  Figure  2a  shows  that  our  wavelet/linear  state  space  method  accurately  reconstructed  the 
signals,  and  precisely  estimated  the  regulatory  strength.  The  results  agree  with  our  expectation  that  Approach 
B  effectively  reduces  the  impact  of  large  amount  of  noise  and  capture  change  patterns,  which  makes  the 
estimation  of  the  regulatory  strength  and  relationship  more  precise. 

This  part  of  the  project  is  still  ongoing;  we  are  testing  Approach  B  on  the  re-profiled  rat  acute  transcriptional 
time  series  of  bolus  administration  of  glucocorticoids. 


Development  of  proposed  dynamic  Bayesian  network  for  genetic  network  reconstruction  (Year  3) 

We  have  also  developed  Approach  A,  dynamic 
Bayesian  network  (DBN),  and  we  are  applying  it  to 
reconstructing  genetic  networks  based  on  the  rat 
muscle  acute  or  chronic  glucocorticoid  transcriptional 
time  series  data  sets.  DBN  is  a  graphical  model  of 
repeated  segments  to  represent  stochastic  process. 

A  DBN  with  hidden  Markov  model  is  illustrated  in 
Figure  3.  DBN  is  the  most  advanced  tool  for 
analyzing  dynamic  process  and  building  dynamic 
probabilistic  model.  We  will  use  DBN  with  the 
Markov  model  to  analyze  the  regulatory  effects  of 
genes  represented  in  the  transcriptional  time  series. 

Note  that  each  node  in  Figure  3  represents  a  vector 
of  variables,  e.g.,  external  stimuli,  hidden  state 
variables,  or  observations. 

DBN  has  a  number  of  advantages.  It  is  suitable  to  model  complex  hierarchical  relationships  and  relationships 
with  nonlinear  properties  using  dynamic  time  series  data.  In  addition,  DBN  does  not  make  any  restrictive 
assumptions  such  as  independency  among  variables  and  linear  interaction.  On  the  contrary,  DBN  faithfully 
evaluates  all  possible  regulatory  network  structures  and  identifies  the  most  optimal  one  that  is  best  supported 
by  the  data.  Therefore,  DBN  does  not  suffer  from  common  problems  existing  in  most  conventional  statistical 
methods,  such  as  confounding  variables. 

Despite  the  important  advantages  of  Bayesian  network  method,  it  has  a  major  challenge,  it  demands  a  large 
sample  size  that  is  often  hard  to  obtain.10 12  To  maintain  a  balance  between  variable  size  and  sample  size  and 
avoid  overfitting,  we  selected  salient  variables  that  have  the  strongest  influence  on  the  outcomes  of  interest 

based  on  the  mutual  information  criterion,  m(X,  Y)  =  £*,yP(x,y) log  [^^-]  ,  where  P(x,y),  P(x),  and  P(y) 

are,  respectively,  the  joint  and  marginal  probabilities  of  variable  X  and  Y,  and  X  and  Y  are  variable  and 
outcome.  If  variable  X  is  independent  from  outcome,  Y,  the  mutual  information  is  zero.  The  larger  the  mutual 
information  is,  the  stronger  the  influence  of  X  on  Y  is.  Because  that  we  are  modeling  a  dynamic  process,  we 

use  the  generalized  mutual  information  m(X,Y )  =  SLiE*  v^t(x<y)l09  \  ttt  J  >  where  Pt(x,y),  Pt(x),  and 

Pt(y)  are,  respectively,  the  joint  and  marginal  probability  of  variable  X  and  Y  at  time  point  t.  This  generalized 
mutual  information  given  equal  weight  to  those  at  every  time  point.  We  are  trying  to  select  a  set  of  variables 
according  to  the  available  number  of  samples  and  time  points  in  the  time  series,  so  that  our  DBN  model  has  a 
sound  statistical  support  by  the  available  samples.  Moreover,  we  have  developed  a  novel  variable  ranking  and 
selection  method,  the  Bayesian  Network  with  Variable  Combination  (BNVC).13  We  will  also  use  BNVC  to  rank 
the  candidate  input  variables  (genes)  and  select  the  variables  that  have  the  strongest  influence  on  the  desired 
outcomes. 


Figure  3.  A  general  diagram  of  dynamic  Bayesian 
network. 
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We  are  using  the  recently  developed  DBN  software  to  analyze  the  rat  muscle  acute  or  chronic  glucocorticoid 
transcriptional  time  series  data  sets.  The  main  results  are  the  transition  matrix  Pt(Xt+1\Xt )  where  Xt  is  the  set 
of  selected  variables  (genes)  at  time  point,  t.  The  transition  matrix  is  the  probabilistic  trajectory,  and  it 
describes  the  regulatory  effect  of  selected  variables  (genes)  at  the  current  time  point  on  the  values  of  variables 
(genes)  at  the  next  time  point.  With  the  transition  matrix,  we  can  derive  genetic  regulatory  networks. 


Progress  of  Specific  Aim  2 

Specific  Aim  2.  Model  the  beneficial  effects  of  daily  bolus  vs.  detrimental  effects  of  chronic  administration  of 
glucocorticoids. 

As  discussed  above,  since  we  recently  re-profiled  the  rat  acute  transcriptional  time  series  set  and  two 
computational  network  reconstruction  algorithms,  we  are  working  on  network  reconstruction  using  the  new  time 
series  and  the  algorithms.  The  results  will  be  needed  to  complete  the  proposed  work  in  Aim  2,  the  work  on  Aim 
2  is  ongoing. 

Identification  of  commonly  shared  canonical  pathways  related  to  DMD  pathogenesis  at  different 
age/disease  stages  across  human  DMD,  and  mouse  and  dog  models  of  DMD  (Year  2) 

Aim  2d  is  to  find  which  genes  potentially  lead  to  asynchronous  regeneration  in  Duchenne  muscle  and  how 
glucocorticoids  rescue  the  issue.  The  more  we  study  on  this  topic,  the  more  we  realize  that  it  is  important  to 
trace  the  upstream  of  the  disease  pathogenesis,  i.e.,  the  genes  that  are  affected  in  the  early  or  presymptomatic 
stage  (such  as  fetus  or  infant)  and  may  become  seeds  leading  to  broader  molecular  deregulation.  Studying 
these  upstream  genes  not  only  will  increase  our  understanding  of  how  the  disease  is  initiated  and  develop,  but 
also  it  is  important  for  design  better  glucocorticoids  regimen. 

Our  lab  has  multiple  muscle  transcriptional  profiling  data  sets  of  human  DMD,  and  animal  DMD  animal  models 
(e.g.,  mouse  mdx,  and  dog  GRMD  (Golden  Retriever  Muscular  dystrophy))  that  contain  subjects  of  different 
age  stages  (see  list  below).  Based  on  these  data,  we  have  identified  a  series  of  commonly  shared  canonical 
pathways  related  to  DMD  pathogenesis  at  different  age/disease  stages  across  human,  mouse  and  dog  (Table 
1).  Two  age  stages  are  matched  across  human,  mouse  and  dog:  presymptomatic/young  stage  (human  0-2 
year-old,  mouse  1  week-old,  dog  4-9  week-old),  postsymptomatic/old  stage  (human  2-10  year-old,  mouse  3 
week-old,  dog  6  month-old). 

List  of  transcriptional  data  sets  used  in  this  analysis: 

1.  Multi-group  human  neuromuscular  diseases:  10  groups,  109  samples,  profiled  on  Affymetrix  Human 
Genome  U133  plus  2.0  microarray.  The  data  set  includes  DMD  and  control  subjects. 

2.  Early  stage  human  DMD  muscle  time  series:  control  and  DMD  fetus  and  infant  subjects,  total  1 1 
samples,  profiled  on  Affymetrix  Human  Genome  U95Av2  microarrays. 

3.  Mdx  mouse  early  muscle  necrosis  time  series:  wild  type  and  mdx  mice,  4  age  groups  (5  days,  2,  3,  and 
6  weeks),  total  24  samples,  profiled  on  Affymetrix  Mouse  Genome  430  v2.0  microarray. 

4.  GRMD  dog  muscle  time  series:  wild  type  and  GRMD  dogs,  3  muscle  types  (vastus  lateralis,  cranial 
sartorius,  long  digital  extensor),  2  age  groups  (4-9  weeks  and  6  months),  total  72  samples,  profiled  on 
Affymetrix  Canine  Genome  2.0  microarray. 


Table  1 .  List  of  commonly  shared  canonical  pathways  related  to  DMD  pathogenesis  at  different  age/disease 
stages  across  human,  mouse  and  dog. 

Production  of  Nitric  Oxide  and  Reactive  Oxygen  Species  GM-CSF  Signaling 
in  Macrophages 

Germ  Cell-Sertoli  Cell  Junction  Signaling  PI3K  Signaling  in  B  Lymphocytes 
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Leukocyte  Extravasation  Signaling 

Actin  Cytoskeleton  Signaling 

LXR/RXR  Activation 

RhoA  Signaling 

Antigen  Presentation  Pathway 

Ephrin  Receptor  Signaling 

Molecular  Mechanisms  of  Cell  Proliferation 

IL-4  Signaling 

Role  of  NFAT  in  Muscle  Hypertrophy 

Cellular  Effects  of  Sildenafil  (Viagra) 

a-Adrenergic  Signaling 

VEGF  Signaling 

Integrin  Signaling 

Complement  System 

Fibrosis 

HGF  Signaling 

P2Y  Purigenic  Receptor  Signaling  Pathway 

FcY  Receptor-mediated  Phagocytosis  in  Macrophages  and 
Monocytes 

Sphingosine-l-phosphate  Signaling 

fMLP  Signaling  in  Neutrophils 

Nitric  Oxide  Signaling  in  Muscle 

IGF-1  Signaling 

Regulation  of  Actin-based  Motility  by  Rho 

All  of  the  listed  pathways  (Table  1)  are  sporadically  or  not  affected  at  the  presymptomatic  stage  in  human, 
mouse  and  dog,  but  they  are  much  more  broadly  affected  at  the  postsymptomatic  stage. 


Development  of  experimental  model  of  repeated  injury  to  induce  asynchronous  regeneration  in  normal 
mouse  muscle  (Year  1, 2  and  3) 
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We  have  made  significant  progress  in 
testing  and  comparing  several  different 
experimental  models  of  multiple  staged 
injections  of  notexin  for  inducing 
asynchronous  (Figure  4).  To  date,  13 
mice  were  bilaterally  injected  in 
gastrocnemii,  with  0,  1,  2,  4,  5,  or  10 
days  between  the  first  and  second 
injections  (at  least  2  mice/group).  The 
mice  were  given  BrdU  on  the  third  day 
after  the  second  injection,  and  were 
sacrificed  13  days  after  the  second 
injection  (the  time  when  regeneration 
due  to  the  second  injection  should 
largely  complete).  The  purpose  of  this 
spacing  was  to  see  if  the  stage  of  the 
regenerating  muscle  from  the  first  injury 
resulted  in  different  effects  on  the 
second  injury.  The  stages  of  the  first 
injury  when  the  second  injection  was 
place  were  myoblast  activation, 

myoblast  proliferation,  myoblast  differentiation,  myotube  fusion,  and  myotube  maturation,  respectively. 
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Figure  4.  Study  design  for  feasibility  test  of  repeated  injections 
for  inducing  asynchronous  regeneration  in  normal  mouse 
muscle.  As  indicated  in  the  figure,  spacing  between  the  two  injections 
was  designed  based  on  interfering  with  specific  stages  of 
regeneration  by  the  second  injection/injury. 


Placement  of  repeated  injection  of  notexin  to  muscles. 

We  learned  a  few  methods  that  are  important  for  giving  proper  repeated  injections  of  notexin  to  muscles.  We 
found  that  dipping  the  needle  in  tattoo  dye  resulted  in  the  surface  of  the  gastrocnemius  being  covered  in  dye. 
The  dye  on  the  muscle  surface  interfered  with  our  ability  to  place  the  successive  injection  accurately  because 
we  could  not  see  the  original  insertion  site  and  only  portions  of  the  injection  tract  through  the  muscle. 
Additionally,  the  tattoo  dye  on  the  surface  of  the  muscle  dried  into  a  layer  of  plastic  that  inhibited  proper  wound 
healing.  Although  we  were  able  to  remove  the  majority  of  this  excess  dye  during  the  second  surgery,  the  skin 


of  some  of  the  animals  did  not  heal  together  properly.  This  may  have  also  caused  an  immune  reaction  and 
slowed  the  muscle  regeneration  -  it  is  not  possible  for  us  to  determine  the  specific  effects  the  excess  tattoo 
dye  had.  As  a  result  of  these  two  complications,  analyzing  the  tissues  has  been  an  extremely  slow  process. 
Importantly,  we  have  found  that  we  cannot  be  sure  if  the  tissue  section  we  are  analyzing  has  been  damaged 
more  than  once.  Based  on  all  of  these  observations,  we  mixed  tattoo  dye  with  the  notexin  and  injected  into  the 
muscle.  While  dye  dispersion  through  the  tissue  may  not  be  to  the  same  extent  as  the  notexin,  the  dye  will 
mark  at  least  a  portion  of  the  tissue  damaged.  Each  round  of  damage  will  be  marked  with  a  different  color 
tattoo  dye  to  enable  us  to  distinguish  regions  of  regeneration.  Including  the  tattoo  dye  in  the  injection  material 
will  guarantee  that  the  area  of  the  tissue  being  analyzed  has  overlapping  regions  of  regeneration  because  we 
will  be  able  to  see  more  than  one  color  of  dye  particles.  Notably,  we  have  tested  the  tattoo  dyes  and  found  that 
they  do  not  interfere  with  fluorescent  immunohistochemical  analysis  of  regeneration. 


Evidences  that  show  repeated  injection  induced  muscle  injuries  and  asynchronous  muscle  regeneration  in 

normal  mouse  muscles 

Using  the  gastrocnemii  isolated  from  the  repeated  injection  area,  we  sectioned  the  tissues  for  histological 
analysis.  We  analyzed  the  images  from  H&E,  Van  Gieson,  and  BrdU  stained  slides,  and  this  work  is  still 
ongoing.  Based  on  the  following  evidences,  our  initial  analysis  showed  that  we  have  successfully  interfered 
with  the  regeneration  process  in  muscles  from  the  groups  where  the  second  injection  occurred  during  a  later 
stage  of  regeneration  (e.g.  10  day  reinjury  group). 

1 .  The  tissues  that  should  be  completely 
regenerated  were  still  undergoing 
rounds  of  regeneration  (Figure  5) 

2.  An  increase  in  degenerating  fibers  was 
found  in  groups  with  more  than  one  day 
between  injections  (Figure  6). 

3.  We  evaluated  fiber  diameter  sizes, 
focusing  solely  on  the  fibers  that  have 
regenerated  at  least  one  time  (we 
measured  only  fibers  with  centralized 
nuclei).  We  found  that  the  10  day 
reinjury  group  had  smaller  fibers  and 
more  of  these  small  fibers  than  any 
other  group,  while  still  maintaining 
larger  fibers.  This  greater  range  in  fiber 
sizes  mimics  what  is  seen  in  patients 
with  DMD  (Figure  6). 

4.  We  used  Van  Gieson  staining  to 
visualize  connective  tissue,  and  to 
determine  if  any  of  the  reinjury  groups 
have  an  increase  incidence  of  fibrosis. 

Cursory  examination  suggests  that  the 


muscle  of  the  reinjury  groups  with  more 
days  in  between  two  injuries  have  an  increase  in  collagen  (Figure  7) 


Figure  5.  Simultaneous  muscle  regeneration  and 
degeneration  observed  in  repeated  injuried  muscles.  This 
gastrocnemius  from  the  day  10  reinjury  model/group  still  had 
regeneration  and  degeneration  at  13  days  after  the  second 
injury.  Normal  muscles  with  single  injury  should  almost  fully 
regenerate  at  this  stage.  However,  in  this  case,  there  is  an 
increase  in  degenerating  fibers  (yellow  arrow)  and  regeneration 
is  still  ongoing:  presence  of  myoblast  fusion  (blue  arrow)  and 
dividing  myoblasts  (green  arrow). 
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5.  Fibrosis  is  a  major  concern  in  DMD  muslces  because  fibroblasts  migrate  into  the  muscle  tissue  and 
replace  the  area  that  had  been  used  to  generate  contractile  force  with  connective  tissue.  In  order  to 
determine  if  the  timing  of  the  reinjury  affects  the  amount  of  fibrosis,  we  performed  Van  Gieson  staining 
to  measure  the  amount  of  interstitial  connective  tissue  present  in  muscle  sections.  We  analyzed  the 
images  of  the  stained  muscle  sections  using  Image  J  software  to  determine  the  percentage  of  muscle 
tissue  per  image  field  that  is  occupied  by  connective  tissue.  We  found  that  the  10  day  reinjury  group 
has  the  greatest  increase  in  connective  tissue  compared  to  other  groups  (Figure  8).  The  5  days  injury 


Figure  7.  Significant  increase  of  fibrosis  due 
to  repeated  injuries.  Considerable  fibrosis  (light 
pink  staining  between  fibers  and  dark  red)  exists 
in  this  image  of  gastrocnemius  from  an  8  wk  old 
BL6  mouse  from  the  day  1 0  reinjury  group. 
Although  the  analysis  of  the  Van  Gieson  data  is 
in  progress,  the  preliminary  data  suggest  there  is 
a  significant  increase  in  fibrosis  in  the  day  10 
reinjury  group  over  groups  with  2  or  fewer  days 
between  injuries. 

group  also  had  an  increase  in  interstitial  connective  tissue.  Next,  we  tried  to  determine  if  this  change  in 
connective  tissue  was  permanent.  When  connective  tissue  becomes  a  permanent  fixture  of  the  muscle 
it  is  considered  fibrosis,  which  is  characterized  by 
a  transition  from  collagen  III,  a  plastic  collagen  that 
is  often  involved  in  remodeling,  to  collagen  I.  To 
determine  if  this  switch  is  occurring,  we  stained  for 
collagen  III  and  collagen  I.  We  found  that  there 
was  some  degree  of  conversion  from  collagen  III  to 
collagen  I  in  all  groups  except  the  0  day  and  10 
day  reinjury,  This  is  understandable  because  the  0 
day  injury  group  regenerates  its  muscle  in  the 
stage-specific  manner  that  produces  successful 
regeneration.  There  is  a  potential  explanation  for 
the  lack  of  collagen  I  in  the  10  day  reinjury  group: 
fresh  rounds  of  the  regeneration-degeneration 
cycle  are  ongoing  -  resulting  in  continued 
remodeling  inhibiting  conversion  of  collagen  III  to 
collagen  I. 

6.  We  have  also  noticed  that  two  of  the  samples  from 

the  10  day  reinjury  group  have  calcified  fibers  (Figure  9).  We  are  planning  to  look  further  into  what  their 
presence  means.  The  calcified  fibers  are  also  seen  in  human  DMD  muscles. 


Figure  8.  There  is  an  increase  in  interstitial 
connective  tissue  in  10  day  reinjury  group,  y-axis 
is  ratio  of  fibrosis  to  total  tissue  area  per  field. 


Figure  6.  Cumulative  frequency  distribution  of 
centrally  nucleated  fibers  shows  that  10  day  reinjury 
group  has  the  smallest  and  largest  range  of  fibers,  y- 

axis  is  cumulative  frequency  distribution,  x-axis  is  Feret’s 
minimum  diameter  (urn). 
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Laser  capture  microdissection  (LCM)  for  sectioning  muscle  regions  with  different  regeneration  stages  induced 

by  repeated  injury 

We  used  LCM  to  section  the  areas  from  staged 
injuries,  and  completed  RNA  and  protein  isolation  from 
the  sectioned  tissues  in  order  to  measure  the  levels  of 
various  markers  of  regeneration,  adipogenesis, 
fibrosis,  and  regeneration.  This  will  further  our 
understanding  of  the  molecules  involved  in  the  cross¬ 
talk  between  neighboring  injuries  and  elucidate  the 
localization  of  these  molecules  within  the  injury 
microenvironments.  Whole  genome  mRNA  expression 
profiling  has  been  performed  on  microdissected 
sections  captured  from  two  injured  adjacent  regions 
and  the  middle  area  to  identify  molecular  signature  in 
each  region  (Figures  10B,  IOC  and  10D)  and  thus 
evaluate  their  involvement  in  cross-talk  between 
neighboring  injury  foci.  We  have  sectioned  a 
representative  of  each  of  the  reinjury  groups  for  LCM. 

These  representatives  were  chosen  because  at  least 
one,  preferably  both,  of  the  tattoo  dyes  that  labeled  the  injection  tracts  was  seen  in  our  aforementioned  H&E 
images.  LCM  was  used  to  isolate  small  areas  of  tissue  associated  with  each  individual  injection,  as  well  as  the 
region  between.  From  these  sections  we  have  purified  good  quality  RNA.  We  are  analyzing  the  mRNA  profiling 
data  determine  if  we  can  elucidate  what  molecules  are  involved  in  crosstalk.  This  work  is  being  partially 
supported  by  other  funding  resources. 


Figure  9.  Calcified  fibers  were  seen  in  10  day 
reinjury  muscle  samples. 


Show  that  daily  glucocorticoids  result  in  more 
successful  regeneration  in  both  normal  asynchronous 

regeneration  and  mdx  mice. 

Based  on  the  above  results  suggesting  that  we  are 
able  to  create  some  of  the  hallmarks  of  DMD  using 
repeated  injury  model,  we  are  planning  to  perform 
Aim  3b  spacing  injuries  using  the  10  day  reinjury 
model.  The  IACUC  protocol  for  this  experiment  has 
been  submitted  and  is  currently  under  review.  We  are 
using  glucocorticoids  in  treating  these  repeated  injury 
groups. 

In  addition,  we  have  conducted  a  study  to  determine 
the  effects  of  time  of  glucocorticoid  administration 
with  respect  to  the  animals’  endogenous  circadian 
rhythm  (we  used  the  light  and  dark  cycle  of  the  animal 
facility  as  our  benchmark).  Animals  were  treated  with 
prednisone  orally  at  lights  ON  (circadian  nadir)  or  lights  OFF  (circadian  peak)  for  five  weeks.  While  we  did  not 
see  a  significant  difference  in  our  functional  assays,  we  found  histologically  that  administering  glucocorticoids 
at  the  animals’  endogenous  circadian  nadir,  rather  than  its  endogenous  peak,  resulted  in  a  significant  increase 
in  fiber  size  and  a  significant  decrease  in  inflammation  and  fibrosis  (Figure  11).  This  result  has  been  included 
in  one  of  our  recent  publications  14 


Figure  10.  Laser  capture  microdissection  of  muscle 
regions  with  different  regeneration  stages. 
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Additional  ongoing 
experiments:  We  will  be 
continuing  our  histological 
characterization  of  these  tissues. 

We  plan  to  assess  the 
inflammation  and  continued 
regeneration-degeneration  seen 
in  the  10  day  reinjury.  We  will  use 
a  general  macrophage  stain  to 
measure  the  degree  of 
inflammation  in  muscle  sections. 

Then  we  will  further  characterize 
the  marchrophages  present  by 
staining  for  markers  of  specific 
macrophage  types  because  some 
macrophages  (type  2)  are 

necessary  for  regeneration.  We  also  plan  to  do  a  stain  specific  for  degenerating  fibers  to  confirm  the  increase 
we  saw  in  the  H&E  counts  we  previously  performed. 


Figure  11.  Effect  of  timing  of  administering  glucocorticoids.  The  muscles 
of  mice  treated  at  the  endogenous  circadian  nadir  (A)  results  in  significantly 
larger  fibers  and  less  inflammation  and  fibrosis  compared  to  treated  at  the 
endogenous  circadian  peak  (B). 


Key  Research  Accomplishments 

1 .  We  have  successfully  developed  and  tested  the  computational  genetic  network  reconstruction 
algorithms  proposed  in  Aim  1.  Using  simulation  data,  we  showed  that  the  methods  are  very  effective  in 
reducing  noise  and  estimating  regulatory  strength  and  relationship  among  genes. 

2.  We  have  re-profiled  rat  time  series  rat  acute  transcriptional  time  series  of  bolus  administration  of 
glucocorticoids  using  lllumina  gene  expression  BeadChip  that  measures  four  times  of  genes  relative  to 
the  old  Affymetrix  chip  used  to  generate  the  existing  time  series.  We  are  using  the  new  lllunimna  time 
series  to  reconstruct  regulatory  networks  as  originally  proposed  in  the  grant. 

3.  By  integrating  our  pre-existing  multiple  muscle  transcriptional  profiling  data  sets  of  human  DMD, 
mouse  mdx,  and  dog  GRMD,  we  have  identified  a  series  of  commonly  shared  canonical  pathways 
related  to  DMD  pathogenesis  at  different  age/disease  stages  across  human,  mouse  and  dog.  All  these 
pathways  (Table  1)  are  sporadically  or  not  affected  at  the  presymptomatic  stage  in  human,  mouse  and 
dog,  but  they  are  much  more  broadly  affected  at  the  postsymptomatic  stage. 

4.  We  have  tested  multiple  injection  models  for  inducing  asynchronous  regeneration  in  normal  mouse 
muscle.  Based  on  the  results  reported  above,  we  found  the  reinjury  model  with  10  days  in  between  two 
injections  of  notexin  lead  to  pathological  damages  that  most  closely  resemble  the  muscle  pathology 
observed  in  mdx.  We  plan  to  use  this  injury  model  to  induce  asynchronous  muscle  regeneration,  and 
then  use  the  injured  muscles  in  testing  the  effect  of  glucocorticoids  in  resynchronizing  regeneration. 

5.  We  have  tested  and  compared  the  effect  of  timing  of  glucocorticoids  administration  with  respect  to  the 
animal’s  endogenous  circadian  rhythm.  Our  preliminary  results  show  significant  histological  difference 
in  terms  of  success  of  muscle  recovery  from  repeated  injuries  in  comparing  muscles  treated  with 
prednisone  orally  at  lights  ON  (circadian  nadir)  or  lights  OFF  (circadian  peak)  for  five  weeks.  This  result 
implies  the  connection  of  the  action  of  glucocorticoid  drug  and  the  circadian  rhythm. 


Reportable  Outcomes 

1.  Evolved  from  this  award,  an  MDA  grant  has  been  awarded,  it  is  a  comprehensive  project  with  emphasis 
on  biological  validation:  Mechanism  of  steroid  action  in  DMD,  grant#  15914,  1/1/2010-  12/31/2012. 
The  goal  is  to  understand  the  effects  of  glucocorticoids  on  muscle  and  enable  the  development  of 
better  targeted  and  more  effective  therapies  for  Duchenne  muscular  dystrophy  dynamically.  This  MDA 
grant  proposal  is  led  by  Dr.  Eric  Hoffman,  and  it  supports  a  doctoral  student's  dissertation  project. 
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2.  The  PI,  Dr.  Zuyi  Wang,  was  awarded  a  grant  in  June  201 1  by  PPMD  (Parent  Project  Muscular 
Dystrophy),  a  private  foundation  dedicating  to  DMD  research.  The  purpose  of  the  grant  is  to  prioritize 
FDA-approved  drugs  for  human  and  animal  testing  for  DMD.  The  work  in  this  grant  has  strong 
connections  with  the  newly  awarded  PPMD  grant  in  many  aspects. 

3.  The  PI,  Dr.  Zuyi  Wang,  was  awarded  a  pilot  grant  in  February  2011  by  Clinical  and  T ranslational 
Science  Institute  at  Children's  National  (CTSI-CN)  funded  by  NIH  CTSA  program.  This  pilot  grant 
focuses  on  studying  glucocorticoids  action  in  treating  asthma.  Chronic  inflammation  is  a  shared  feature 
in  many  disorders  of  various  tissues  (including  muscle  and  lung)  for  which  glucocorticoids  provide  the 
standard  of  care.  The  Pi’s  research  on  glucocorticoid  action  in  muscle  has  provided  valuable  resources 
of  experiences  and  methods  to  the  proposed  study  on  asthma. 

4.  The  PI,  Dr.  Zuyi  Wang,  received  two  grants  from  the  Junior  Faculty  Fund  program  of  Center  for  Genetic 
Medicine/Department  of  Integrative  Systems  Biology  in  201 1  and  2012  respectively.  The  two  grants  are 
designed  to  establish  in  vitro  model  for  studying  the  effects  of  glucocorticoid  drugs  in  human  asthmatic 
airway  epithelial  cells.  This  is  a  systems  biology  study  of  molecular  effects  on  glucocorticoid  drug  action 
in  asthma. 

5.  L.  Fu,  R.  Freichtat,  S.  Teach,  R.  Lorenzo,  E.  Hoffman,  Z.  Wang.  “Bayesian  network  analysis  suggests 
ranking  of  major  clinical,  behavioral  and  environmental  variables  influencing  asthma  outcomes.” 
American  Thoracic  Society;  2012;  San  Francisco,  Am  J  Respir  Crit  Care  Med ;  2012. 

6.  A.  Baudy,  E.  Reeves,  J.  Damsker,  C.  Heier,  L.  Garvin,  J.  McCall,  S.  Rayavarapu,  Z.  Wang,  J. 
VanderMeulen,  A.  Sali,  V.  Jahnke,  S.  Dugue,  D.  DuBois,  M.  Rose,  K.  Nagaraju,  and  E.  P.  Hoffman.  “A- 
9,1 1  steroidal  compounds  dissociate  NF-kB  transrepression  from  GRE-mediated  transactivation.”  J. 
Pharmacology  and  Experimental  Therapeutics,  2012,  in  press. 


Conclusion 

We  developed  the  two  proposed  computational  network  reconstruction  methods,  linear  state  space  model  and 
DBN,  to  infer  transcriptional  networks  using  the  rat  acute  transcriptional  time  series  of  bolus  administration  of 
glucocorticoids.  We  have  re-profiled  the  samples  of  the  rat  acute  transcriptional  time  series  of  bolus 
administration  of  glucocorticoids  using  lllumina  gene  expression  BeadChip  to  obtain  a  much  larger  number  of 
measured  genes.  It  is  well  recognized  that  reconstructing  large  scale  genetic  network  is  a  very  complex  task. 
One  major  obstacle  is  that  the  transcriptional  microarray  data  often  contain  considerable  noise  that  may  lead  to 
false  estimation  of  regulatory  relationships  among  genes.  We  are  testing  and  comparing  the  two  developed 
methods  on  the  transcriptional  time  series  data  sets  about  the  performance  in  reliably  and  correctly 
reconstructing  the  regulatory  networks.  This  work  is  ongoing. 

Through  integrating  multiple  muscle  data  sets  of  human  DMD,  mdx  and  GRMD,  we  have  identified  a  series  of 
commonly  shared  canonical  pathways  related  to  DMD  pathogenesis  at  different  age/disease  stages  across 
human  DMD,  mouse  and  dog  models  of  DMD.  These  canonical  pathways  are  very  differentially  affected  at  the 
presymptomatic  and  postsymptomatic  stages,  we  are  studying  these  pathways  and  relevant  genes  to  gain  in 
depth  understanding  of  disease  progression. 

We  have  tested  multiple  injection  models  for  inducing  asynchronous  regeneration  in  normal  mouse  muscle. 
Based  on  the  analysis  of  injured  muscles,  we  found  the  reinjury  model  with  10  days  in  between  two  injuries 
lead  to  pathological  damages  that  most  closely  resemble  the  muscle  pathology  observed  in  untreated  mdx 
muscles  and  human  DMD  muscles.  We  are  using  the  10  day  reinjury  model  to  introduce  asynchronous 
regeneration  in  normal  mouse  muscle  and  study  the  effect  of  glucocorticoids  in  resynchronizing  the 
regeneration. 
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Finally,  based  on  part  of  the  work  accomplished  in  this  grant,  we  have  been  awarded  several  new  grants 
focusing  on  studying  glucocorticoid  mechanism  in  treating  DMD  and  asthma,  and  published  one  journal  article 
and  one  conference  abstract. 
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